Multi-dimensional stability of planar Lax shocks in 
hyperbolic-elliptic coupled systems 



Toan Nguyen-^ 



Abstract 

We study nonlinear time-asymptotic stability of small-amplitude planar Lax shocks 
in a model consisting of a system of multi-dimensional conservation laws coupled 
with an elliptic system. Such a model can be found in context of dynamics of a gas 
in presence of radiation. Our main result asserts that the standard uniform Evans 
stability condition implies nonlinear stability. The main analysis is based on the 
earlier developments by Zumbrun for multi-dimensional viscous shock waves and by 
Lattanzio-Mascia-Nguyen-Plaza-Zumbrun for one-dimensional radiative shock pro- 
files. 



1. Introduction 

In the present paper, we consider the following general hyperbolic-elliptic coupled 
system, 

«t + Ei=i/j(^)^, +^divg = 0, ^ 
- V div g + g + Vg{u) = 0, ^ ' ^ 

consisting of a system of conservation laws coupled with or regularized by an elliptic 
system, with imposed initial data u{x,0) = Uo{x). Here, x G M*^, L is a constant 
vector in M", the unknowns m G M" and q G M*^, for n > 1, d > 2, the nonlinear 
vector- valued flux fj{u) G M", and the scalar function g{u) G M. 

The study of (II. ip is motivated by a physical model or a so-called radiating gas 
model that describes dynamics of a gas in presence of radiation. Such a model (due to 
high-temperature effects) consists of the compressible Euler equations coupled with 
an elliptic system representing the radiative flux. See, for example, [s], [26|, for its 
derivations and discussions further on physical applications. 

The system (11.11) in its spatially one-dimensional form has been extensively stud- 



ied by many authors such as Liu, Schochet, and Tadmor [23|, [ISJ, Kawashima and 
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Nishibata 0, M, EoLSerre [25|, Ito 0, Lin, Coulombel, and Goudon 0, [l6 
among others. In Lattanzio, Mascia, and Serre show the existence and regular- 
ity of (planar) shock profiles (whose precise definition will be recalled shortly below) 
in a general setting as in fll.lh . and recently in a collaboration with Lattanzio, Mascia, 
Plaza, and Zumbrun IJ, |2l|, we show that such radiative shocks with small ampli- 
tudes are nonlinearly asymptotically orbitally stable. Regarding asymptotic stability, 
all of aforementioned references deal with spatially one-dimensional perturbations. In 
this work, we are interested in asymptotic stability of such a shock profile with respect 
to multi-dimensional perturbations. Regarding asymptotic behaviors of solutions to 
the model system (11. ip in the multi-dimensional spaces, we mention recent related 



works by Wang and Wang [27[ and by Liu and Kawashima [17[. There, however, the 
authors study stability of constant states (or the zero state) and the model system 
(II. ip that they consider is restricted to the case when u are scalar functions. In this 
paper, we study stability of planar shocks and allow u to be vector-valued functions. 

1.1. Shock profiles 

To state precisely the objective of our study, let us consider the one-dimensional 
system of conservation laws: 

ut + fiiu),,=0, (1.2) 

for vector function u G M". We assume that the system is strictly hyperbolic, that 
is, the the Jacobian matrix dfi{u) has n distinct real eigenvalues Xj{u), j = 1, ■ ■ ■ ,n, 
with Ai(m) < ■ ■ ■ < Xn{,u), for all u. It is easy to see that such a conservation laws 
(II. 2p admits weak solutions of the form u = u{x — st) with 



M+, X > Xo, 

U-, X < Xq, 



for u± G M"', s G M, and Xq G R, assuming that the triple {u±, s) satisfies the 
Rankine-Hugoniot jump condition: 

/i(m+) -/i(m-) = s(m+-m_). (1.3) 

Here, by translation invariant, we take xq = 0. The triple {u±,s) is then called 
a hyperbolic shock solution of the system (ll.2p . It is called a hyperbolic p- 
Lax shock solution of (11.20 if the triple further satisfies the classical p-Lax entropy 
conditions: 

Ap(n+) < s < Ap+i(M+), ^^^^ 
Ap_i(n_) < s < Ap(n_), 

for some p such that 1 < p < n. 
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Next, let us consider the one-dimensional hyperbolic-elliptic system, that is the 
system (11.21) coupled with an elliptic equation: 



-qU + + 9{nU = 0, ^'-"^ 

for vector functions u G M" and scalar G M. Lattanzio, Mascia, and Serre have 
shown ([lil) that there exist traveling wave solutions of (11 .Sp that associate with (or 
regularize) the hyperbolic p-Lax shock. To recall their result more precisely, let us 
denote Lp{u) and Rp{u) the eigenvectors of dfiiu) associated to the eigenvalue \p{u). 
Assume also that the p^^ characteristic field is genuinely nonlinear, that is, 

{VXpY -Rp^O, (1.6) 

and furthermore at the end states u±, there holds the positive diffusion condition 

Lp{u±){Ldg{u±))Rp{u±) > 0. (1.7) 

Here, dg{u±) is the Jacobian row vector in R", consisting the partial derivatives in 
Uj of g{u). The condition (I1.7P indeed comes naturally from the Chapman- Enskog 



expansion, giving a right sign of the diffusion term; see, for example, [23| or [12 



We recall the result in 12 



Given a hyperbolic p-Lax shock {u±, s) of ( 11. 2p and the assumptions ( 11. 6p and 
( ll.7p . there exists a traveling wave solution {u, q^) of (11.51) with the same speed s and 
with asymptotic constant states {u±,0): 

iu,q'){xi,t) = iU,Q'){xi-st), (f/,gi)(±cx)) = (m±,0), (1.8) 

Furthermore, when the shock has a sufficiently small amplitude: — u_\ ^ 1, the 
traveling wave solution is unique (up to a translation shift) and regular (see Theorems 



1.6 and 1.7 of fl2] for precise and much more general statements). 



We call such a traveling wave (II. 80 a radiative p-Lax shock profile. Let Q = 
{Q\0) G M.'^. It is clear that ([/, Q) is a particular solution to the multi-dimensional 
hyperbolic-elliptic system (II. ip . with {U,Q^) as in (ll.8p . We then call the solution 
([/, Q) the planar radiative p-Lax shock of (II. ip . Without loss of generality (that 
is, by re-defining /i by /i — su), in what follows we assume that the shock speed s is 
zero. 

In this paper, we study nonlinear time- asymptotic stability of such a planar radia- 
tive p-Lax shock {U,Q) with sufficiently small amplitudes: \u+ — m_| <^ 1. We shall 
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make several technical and structural assumptions. Our first set of assumptions, as 
a summary of the above assumptions, reads as follows: 

(51) The system (ll.2p is strictly hyperbolic, and the triple {u±, 0) is a hyperbolic 
p-Lax shock of ( II. 2p . 

(52) The system (11. 5p satisfies the genuine nonlinearity and the positive diffusion 
conditions fll.6p and fll.7p 

By hyperbolicity, it is straightforward to see that as long as the shock profile 
{U, Q) is smooth, it enjoys the exponential convergence to their end states, precisely. 



as — )■ +00, for some > 0, > 0. See, for example, a simple proof in |14| |. 
Lemma 2.1. 

In addition, we remark that the condition (11.41) implies that \p{U {xi)) must vanish 
at some point G M along the shock profile. By translating Xi to Xi + x\, we 
assume that it vanishes at Xi = 0. We call such a point singular simply because the 
associated ODE system obtained from the standard resolvent equations is singular at 
this point. For further discussions on this point, see the paragraph nearby equation 
f ll.l7p . Throughout the paper, we assume that 

(S3) xi = is the unique singular point such that Ap(f/(0)) = 0. Furthermore, 
at this point, we assume 



The uniqueness assumption is purely for sake of simplicity. The case of finite 
numbers of singular points should follow similarly from our analysis. 

1.2. Structural assumptions 

We shall make our second set of assumptions on structure of the system (11. ip . Let 
us recall that dfj and dg denote the Jacobians of the nonlinear flux functions fj and 
g, respectively. Let U be some neighborhood in R" of the shock profile U, constructed 
in the previous subsection. Our next assumption concerns the symmetrizability of 
the system. 

(Al) There exists a symmetric, positive definite Aq = Aq{u) such that AQ{u)dfj{u) 
is symmetric and AQ{u)Ldg{u) is positive semi-definite, for all u ^U. 



{d/dxi)\U -u±,Q) < Ce""! 



(1.9) 



— A,(C/(xO)u,.o^O. 



(1.10) 
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One may notice that (Al) is a common assumption in the stabihty theory of 
conservation laws, which may go back to the original idea of Godunov and Friedrichs 
(see, e.g., ^). Essentially, by the standard symmetrizer or W energy estimates, 
the assumption (Al) yields the necessary local well-posedness, and is closely related 
to existence of an associated convex entropy of the hyperbolic system. 

We next impose the well-known Kawashima and Shizuta (KS) condition, which 
has played a very crucial role in studies of time-asymptotic stability. The assumption 
reads 

(A2) For each ^ {0}, no eigenvector of C,jdfj{u±) lies in the kernel of 

Our use of the (KS) condition is to derive sufficient H^, for large s, energy esti- 
mates, and therefore provide sufficient control of "high-frequency" part of the solution 
operator. Here and in what follows, by high- or low-frequency regions, we always mean 
the regions at the level of resolvent solutions that |(A,^)| is large or small, with (A,^) 
being the Laplace and Fourier transformed variables of time t and the spatial variable 
X transversal to Xi. 

1.3. Technical hypotheses at hyperbolic level 

Along with the above structural assumptions, we shall further make the following 
two technical hypotheses at the hyperbolic level (i.e., the level without the presence 
of q in our model (11. ip ). 

(HI) The eigenvalues of J2j ^jdfj{u±) have constant multiplicity with respect to 
^ G ^ ^ 0. 

(H2) The set of branch points of the eigenvalues of {dfi)^^{ir + '^^jdfj){u±), 
r e M, ^ G M°'~^ is the (possibly intersecting) union of finitely many smooth curves r = 
?7^(^), on which the branching eigenvalue has constant multiplicity Sq (by definition 
>2). 

These hypotheses are crucially used in our construction of the Green kernel and 
the resolvent solution in the low-frequency regimes, and borrowed directly from the 



previous analyses introduced by Zumbrun |28|, |29|. The condition (HI) is the stan- 
dard non-strict hyperbolicity with constant multiplicity assumption. Whereas, (H2) 
concerns singularities of the branching eigenvalues. It provides certain compactness 
properties that allow us to later on perform matrix perturbations with acceptable 
errors. We refer the interested reader to ^8|, Section 4.1, for a thorough discussion 
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on these conditions. In particular, (H2) is satisfied always in dimension d = 2 oi for 
rotationally invariant systems in dimensions d > 2. 

It is perhaps worthwhile to mention that these hypotheses might be weakened 
or dropped as observed in [20] for the case of hyperbolic-parabolic settings. More 
precisely, we were able to allow eigenvalues with variable multiplicities (for instance, 
in case of the compressible magnetohydrodynamics equations) and to drop or remove 
the technical condition (H2) in establishing the stability. However, we leave it for 
the future work, as our current purpose is to show that the well-developed stability 



theory [28|, l29[ for the hyperbolic-parabolic systems can be adapted into the current 
hyperbohc-elliptic settings despite the presence of singularity in the eigenvalue ODE 
systems, among other technicalities. 

Finally, regarding regularity of the system, we make the following additional as- 
sumption: 

(HO) fj,g,Ao e C"+\ for some s large, s > s{d) with s{d) := [{d - l)/2] + 5. 

The regularity is not optimal due to repeated use of Sobolev embeddings in our 
estimates of the solution operator, especially the energy-type estimate of the high- 
frequency solution operator in Section H31 One could lower the required regularity by 
deriving much more detailed description of the resolvent solution following Zumbrun 



29|, instead of using the energy-type estimate, in the high-frequency regime. 



Throughout the paper, to avoid repetition let us say Assumption (S) to mean the 
set of Assumptions (SI), S(2), and (S3); Assumption (A) for (Al) and (A2); and. 
Assumption (H) for (HO), (HI), and (H2). 

1.4- The Evans function condition 

As briefly mentioned in the Abstract of the paper, we prove a theorem assert- 
ing that an Evans function condition implies nonlinear time- asymptotic stability of 
small radiative shock profiles , under Assumptions (S), (A), and (H) mentioned earlier. 
Shortly below, we shall introduce the Evans function condition that is sufficient for 
the stability. To do so, let us formally write the system (11. ip in a nonlocal form: 

f ut + Yfj=i fj{u)xj - LdiYlCVg{u) = 0, 

with /C := (— Vdiv ■ +1)~^. We then linearize the system around the shock profile U. 
The linearization formally reads 

Ut-Cu = 0, Cu:= -Y.AAj{xi)u)^^- Ju (1.12) 
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with initial data u{0) = uo, and J'u := —Ldw}C'V{B{xi)u). Here, we denote 
Aj{xi) := dfj{U{xi)) and B{xi) := dg{U{xi)). Hence, the Laplace-Fourier trans- 
form, with respect to variables x the transversal variable, applied to equation 
ffLT2D gives 

\u-C^u = S (1.13) 

where source 5* is the initial data uq. An evident necessary condition for stability is 
the absence of solutions for values of A in {3fJeA > 0}, for each ^ G M^"^, noting 
that, when ^ = 0, A = is the eigenvalue associated to translation invariance. 

We establish a sufficient condition for stability, namely, the strong spectral sta- 
bility condition, expressing in term of the Evans function. For a precise statement, 
let us denote -D±(A,^) (see their definition in f l2.32p below) the two Evans functions 
associated with the linearized operator about the profile in regions Xi ^ 0, corre- 
spondingly. Let C = (f. A). Introduce polar coordinates C, = pC, with ( = (f , A) on 
the sphere S'^, and write D±{X,C^) as D±{(,p). Let us define S'^ = H {3fteA > 0}. 
Our strong spectral (or uniform Evans) stability assumption then reads 

(D) D±{(,p) vanishes to precisely the first order at p = for all C G 5^ and has 
no other zeros in 5*1 x M , . 



The assumption is assumed as in the general framework of Zumbrun |28l . |29 . 
Possibly, it can be verified for small-amplitudes shocks by the work of Freistiihler and 
Szmolyan jif . It is also worth mentioning an interesting work of Plaza and Zumbrun 



22[, verifying the assumption in one-dimensional case. In addition, the assumption 
can also be efficiently numerically checkable; see, for example, numerical computations 
in joj for the case of gas dynamics. 

We remark that even though we only consider in this paper the strong form 
of the spectral stability assumption (D), in the same vein of the main analysis in 
(isl . [iol , our results should hold for a weaker form (thus more precise description for 
stability), namely, the refined stability assumption which involves signs of the second 
derivatives of D±{(, p) in p. In addition, extensions to nonclassical shocks should also 
be possible. Nevertheless, we shall omit to carry out all these possible extensions 
and confine the presentation to the case of the classical Lax shocks under the strong 
spectral assumption (D). 

1.5. Main result 

We are now ready to state our main result. 

Theorem 1.1. Let {U, Q) be the Lax radiative shock profile. Assume all Assumptions 
(S), (A), (H), and the strong spectral stability assumption (D). Then, the profile {U, Q) 
with small amplitude is time- asymptotically nonlinearly stable in dimensions d >2. 
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More precisely, let {u, q) he the solution to (II. ip with initial data uq such that 
the initial perturbation Uq := Uq — U is sufficiently small in L} fl , for some s > 
[{d — l)/2] + 5. Then {u,q){t) exists globally in time and satisfies 



d-l , , 



\u{x, t) - f/(xi)Up < C{1 + t)-^(i-iM+e 
\u{x,t) - U{x,)\hs < C{1 + t)-'^''-'y'\uo\L^nHs 



and 

\qix,t) - Qix,)\w^,, < Cil + ty'^^'-'M+'luolL^r^Hs 
\q{x,t) - Q{xi)\hs+i < C{1 + t)-'/^\uo\LinH' 
for all p > 2; here, e > is arbitrarily small in case of d = 2, and e = when d > 3. 

We obtain the same rate of decay in time as in the case of hyperbohc-parabohc 
setting (see, e.g., (i^). This is indeed due to the fact that in low-frequency regimes 
the estimates for the Green kernel for both cases, here for the radiative systems and 
there for the hyperbolic-parabolic systems, are essentially the same, away from the 
singular point occurring in the first-order ODE system for the former case. 

Let us briefly mention the abstract framework to obtain the main theorem. First, 
we look at the perturbation equations with respect to perturbation variable u = u — U, 
namely, 

ut - Cu = N(u,u^)^, (1.14) 

where Cu = — ^j{Aj{xi)u)xj — Ju as defined in (11.120 and N{u, Ux) is the nonlinear 
remainder term. Since Cu is a zero-order perturbation of the generator — 'Yliji^ji^i)'^)^ 
of a hyperbolic equation, it generates a C° semigroup e^* on the usual space which 
enjoys the inverse Laplace-Fourier transform formulae 



1 /'7H-200 n 

/(^) = ^r-W / / - C^)-'f{xiA)d^d\ (1.15) 

[2mY J^_i^ Jjjd-i 

where C^ is the Fourier-transformed version of the operator C in the transversal 
variable x. 

Having the solution operator e"^* expressed as in (ll.lSp . we may now write the 
solution of (I1.14P by using Duhamel's principle as 

u{x,t)=e'^'uo{x)+ [ e'^'''-'^N{u,Ux)x{x,s)ds, (1.16) 
Jo 

noting that q can always be recovered from u by q{x,t) = —)C(yg{u)){x,t). Hence, 
the nonlinear problem is reduced to study the solution operator at the linearized level, 
or more precisely, to study the resolvent solution of the resolvent equation 

(A - C,)u = f. 
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The procedure might be greatly comphcated by the circumstance that the resulting 
(n + 2) X (ri + 2) first-order ODE system 



(e(xi)W^)., = A{x„X,OW, e(xi) := f^^lj"^^) ° ) , (1.17) 



is singular at the point where the determinant of A^{xi) vanishes, with B dropping 



from rank n + 2 to n + 1. However, as already observed in |14j . |21[ |. we find in the 



end as usual that the Green kernel ^^j^^^ can be constructed, and contribution of the 
terms due to the singular point turns out to be time-exponentially decaying. 

The paper is organized as follows. In Section 2, we will study the resolvent so- 
lutions in low-frequency regions and define the two Evans functions, essential to the 
derivation of the pointwise Green kernel bounds which will be presented in Section 
3. Once the resolvent bounds are obtained, estimates for the solution operator are 
straightforward, which will be sketched in Section 4. A damping nonlinear energy 
estimate is needed for nonlinear stability argument, and is derived in Section 5. In 
the final section, we recall the standard nonlinear argument where we use all previous 
linearized information to obtain the main theorem. 

2. Resolvent solutions and the two Evans functions 

In this section, we shall construct resolvent solutions and introduce the two Evans 
functions that are crucial to our later analysis of constructing the resolvent kernel. 
We consider the linearization of (11. ip around the shock profile {U, Q) 

d 

ut + ^{Aj{xi)u)^^ + Ldiv g = 0, 
i=i 

-V div g + g + V(5(xi)m) = 0, 

where Aj{xi) = dfj{U{xi)), B{xi) = dg{U{xi)), and q = {q^,q'^, ■ ■ ■ ,q'^) G M'^. Since 
the coefficients depend only xi (through U{xi)), we can apply the Laplace-Fourier 
transform to the system (12.11) in time t and transversal variables x. Let us ignore for 
a moment the contribution from the initial data. The Laplace-Fourier transformed 
system then reads 

(A + iA^)u + (Aiu)^, + Lql^ + iLq^ = 0, 

-{qi+tq%,+q' + iBu),,=0, (2.2) 
-iCj {ql, +iq^) + q' + tijBu = 0, ] ^ I, 
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where for simplicity we have denoted := 1^ '■— Multi- 

plying the last equations by i^j, j 7^ 1, and summing up the result, we obtain 

iql,+^q^)\i\' + ^q^-\i\'Bu = 0. 

From this identity, we can solve iq^ in term of u and and then substitute it into 
the first two identities in the system (12.21) . We then obtain 

V 1 + |^|2 J ^ ' 1+|^|2 (2.3) 

-qI,., + (1 + ieT)g' + (Bu),. = 0. 

System (12. 3p is a simplified and explicit version of our previous abstract form Xu — 
C^u = 0, where £| is defined as the Fourier transform of the linearized operator C. 
Now, by defining 

pi := Bu-ql^, 

we then easily derive the following first order ODE system from (12. 3p 

(A.u)^, = -(A + lA^ + LB)u + (1 + leT)"' V, 

ql^=Bu-p\ (2.4) 

= -(1 + 

The key observation here is that this first-order ODE system is very similar to the 
system that we have studied for the one-dimensional case, considering the variable ^ 
parameter. 

2.1. Stable/unstable dimensions 

Next, we can diagonalize Ai with recalling that Ai{xi) = dfi{U{xi)) has distinct 
and nonzero eigenvalues by hyperbolicity assumption (SI). Let us denote ap{xi) = 
Xp{U{xi)) with Xp{U) being the p^^ eigenvalue of dfi{U), introduced in Section [TTTl 
By hyperbolicity, there exists a bounded diagonalization matrix T{xi) such that the 
matrix Ai{xi) can be diagonalized as follows: 



ii(xi) :=r-%T(xi) = ap{xi) (2.5) 




where a_ is the {p — 1) x (p— 1) matrix and negative definite, a+ is the (n—p) x [n—p) 
matrix and positive definite, and Op G M, satisfying ap(+oo) < < ap(— 00) (by the 
Lax entropy conditions (11.41) ). 



10 



Defining v := T ^u, we thus obtain a diagonalized system from (12.4p : 

(Aiv),, = - (a + tA^ + LB + {T-'),,A,Ty, + (1 + \^\yLp\ 

ql^=Bv-p\ (2.6) 

pI = -(1 + iir)g\ 




wliere L := T'^L, B := BT, and = T'^A^T. 

We shall construct the Green kernel for this diagonalized ODE system (12. 6p . To 
do so, let us write tlie system fl2.6p in our usual matrix form with unknown W := 

{Q{x,)W),^^=A{x,,X,OW, (2.7) 

where we have denoted 0(xi) := j ^ J and 

/-{\ + iA^ + LB + {T-^),,AiT) 
A(xi,A,0 := B 

V -(1 + m 

We note that since ap(0) = (see the assumption (S3)), the matrix Ai, and thus 9, 
is degenerate at xi = 0. We shall see shortly below that this singular point causes 
the inconsistency in dimensions of unstable and stable manifolds, and thus the usual 
definition of the Evans function must be modified. 
Let us denote the limits of the coefficients as 

i± := lim A{xi), B± := lim B{xi), L± := lim L(xi), (2.8) 

and 

( -A^\X + tAi + L±B±) (l + |fP)-M±iL±\ 

A±(A,|):= B± -1 . (2.9) 

V -(1 + ieT) ; 

Here, note that A±(A,^) is not quite the limiting matrix of A(xi,A,.^) at infinities. 
Having defined these asymptotic matrices, the asymptotic system of (12. 7p can be 
written as 

W = A±{X,i)W. (2.10) 

We need to determine the dimensions of the stable /unstable eigenspaces. Observe 
that simple computations show 

det(/i - A±) =/i2 det(/i + i±^(A + iAi + L±B±)) - (1 + |f H det(;U + i±^(A + iAi)), 
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where since for p = |(A,f)| — )■ the absolute value of ^±^(A + iA^.) = 0{p), the 
above yields one strictly positive and one strictly negative eigenvalues at each side 
of a; = ±00, denoting pf and ^^+2 (later on, giving one fast-decaying and one fast- 
growing modes). Looking at slow eigenvalues p = 0{p), one easily obtains that 
the first term in the above computation of det(/x — A±) contributes 0{p^) and thus 
eigenvalues p, are of the form 

pf{\l)=P%{\l) + 0{p'), (2.11) 

where /i^ are eigenvalues of —A^{X + iAj_). Now, notice that A^{\ + iAj_) has no 
center subspace, i.e., no purely imaginary eigenvalue, for 3fJeA > 0. Indeed, if it were 
one, say 2^1, then A^{\ + iA^j^)v = i^iv, or equivalently, \v = — X]y=i ^O^i"^' 



some V G M", which shows that A G by hyperbolicity of the matrix Yl'i=i^j^±- 



Thus, A^^{X + iA^.) has no center subspace. Consequently, the numbers of sta- 
ble/unstable eigenvalues of ^^^(A + iA^.) persist as |f| — )■ 0, and remain the same 
as those of A'^^, and thus of A}^. We readily conclude that at x = +cx3, there are 
p + 1 unstable eigenvalues (i.e., those with positive real parts) and n — p + 1 stable 
eigenvalues (i.e., those with negative real parts) . The stable S~^{X,^) and unstable 
[/+(A, ^) manifolds, which consist of solutions that decay or grow at +00, respectively, 
have dimensions 

dimf/+(A,0 =P+ 1, 

^ J (2.12) 

dim S'^{X,C,) = n — p + 1, 

in 9ftA > 0. Likewise, there exist n—p+2 unstable eigenvalues and p stable eigenvalues 
so that the stable (solutions which grow at —00) and unstable (solutions which decay 
at —00) manifolds S~{X,C,) and f/~(A,^), respectively, have dimensions 

dimf/~(A,f) = p, 

_^ ^ (2.13) 

dimS" {X,^) = n — p + 2. 

Remark 2.1. Notice that, unlike customary situations in the Evans function liter- 
ature (see, e.g., [H, 28, |2^), the dimensions of the stable (resp. unstable) manifolds 



5*^ and S (resp. and U ) do not agree. 

2.2. Asymptotic behavior 

We study the asymptotic behavior of solutions to the first order ODE system f l2.7p 
away from the singularity point a: = 0. To simplify our presentation, we consider the 
case when x — )■ +00. Note that our treatment will be unchanged if there were finitely 
many singular points. We pay special attention to the small frequency regime, p — )■ 0. 
By performing a column permutation of the last two columns in (12.91) . with an error of 
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order 0{p'^), and by further performing row reductions with observing that spectrums 
of the two matrices A^^{\ + iA\_) are of order 0{p), strictly separated from ±(9(1), 
we find that there exists a smooth matrix V{\,^) such that 

V-'KV=(^, 2,1 (2.14) 



^0 

with blocks P = diag{P+, P_} + 0{p) with ±3fJeP± > and 

H,iX,0 ■■= -A-\X + tAi) = -T-'A-\X + tAi)T, 

for T being the diagonalization matrix defined as in f l2.5p . We note that H which 
determines all slow modes is spectrally equivalent to 

-A-\X + tAi) + 0{p^). 
We then obtain the following lemma. 

Lemma 2.2. For p sufficiently small, the spectral system f l2.10p associated to the 
limiting, constant coefficients asymptotic behavior of (12.41) has a basis of solutions 

e4(^^'^)^^V^{^X,i), x^O, j = l,...,n + 2, 

where {V^}, necessarily eigenvectors of A±, consist of 2n slow modes associated to 
slow eigenvalues (as in (12.111) ) 

pf{X,0 = /ifo(A,0 + 0{p') J = 2, ...,n + 1, (2.15) 
with /i^ eigenvalues of —A^^{X + iA^.), and four fast modes, 

pf{x,o = ±Oi+o{p), 

/in+2(A,0 = T^^+2 + (^(p)- 

where 9f and 0^^2 ^'^^ positive constants. 

Proof. As discussed above, there are one eigenvalue with a strictly positive real part 
and one with a strictly negative real part at a; = ±oo, giving four fast modes. 
Whereas, 2n slow modes are determined by the matrix H, which is spectrally equiv- 
alent to 

-A^\X + tAi) + Oip'), 
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which gives the expansion fl2.15p . Constructing the eigenvectors Vp of A± associ- 
ated to these slow eigenvalues can be done similarly as in ^9|, Lemma 4.8, since 

the governing matrix —A~^{\ + iJ^±) is precisely the same as those studied in the 
hyperbolic-parabolic systems. Note that these matrices purely come from the hyper- 
bolic part of the system. 

The main idea of the construction is to use the assumption (HI) to separate the 
slow modes into intermediate-slow (or so-called elliptic) modes for which |3fJe/i^| ~ p, 
super-slow (hyperbolic) modes for which |3fJe/i^| ~ and 53mA is bounded away 
from any associated branch singularities ^7^(0) and super-slow (glancing) modes for 
which |3fJe/i^| ~ (? and is within a small neighborhood of an associated branch 
singularity Tij[^. Finally, thanks to the assumption (H2), the glancing blocks can 
also be diagonalized continuously in A and ^, and thus associated eigenvectors can be 



constructed. We refer to 28|, Lemma 4.19] for details. □ 



In view of the structure of the asymptotic systems, we are able to conclude that 
for each initial condition > 0, the solutions to ( 12. 4p in x\ > xq are spanned by 
decaying/growing modes 

as Xi — > +00, whereas for each initial condition xq < 0, the solutions to (12. 4p are 
spanned in xi < —Xq by growing/decaying modes 

: = {^r,...,^/^-_p+2}, ^^^^^ 

: = {C-p+3'-->0n+2}> 

as Xi — oo. Later on, these modes will be extended on the whole line Xi G M, by 
writing them as linear combinations of the corresponding modes that form a basis of 
solutions in respective regions x < —xq, x > xn, or \x\ < \xq\. 



We rely on the conjugation lemma of [19|| to link such modes to those of the 
limiting constant coefficient system fl2.10p . 

Lemma 2.3. Lemma 4-19] For p sufficiently small, there exist unstable/stable 

(i.e., growing/decaying at +oo and decaying/growing at — oo ) solutions 'ijj^{xi, A, ^), (f)f{xi, A, ^), 

in xi ^ ±xo, of class in xi and continuous in A,^, satisfying 



^f{x,, A,0 = 72i,^±(A,|)e^"(^'«>^V^/(A,0(l + O(e-^l-^l)), 
0f (xi, A,0 = 72i,,±(A,0e^?(''«^^K^(A,e)(l + 0{e-^^^^^)), 



(2.18) 
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where i] > is the decay rate of the traveling wave, and /x^ and Vp are as in Lemma 
\2.S\ above. Here, the factors 

721,^±,721,0± ~ 1 

for fast and intermediate- slow modes, and for hyperbolic super-slow modes, and 

72i,^± ~ 1 + [p-'\'^mX - r/f (01 + p]V 

72i,<^± ~ 1 + [p-'\'^mX - r/f (01 + pf'f 

for glancing super-slow modes, for some t^± , t^± < 1 depending on sj . Here, the 
symbol ~ means that we can obtain upper and lower bounds independent of smallness 
of P- 

Remark 2.4. The factors 791 ^± are viewed as diagonalization errors which were 

introduced by Zumbrun in his study of shock waves for hyperbohc/parabohc systems; 
see Lemma 4.19, 28|, or Lemma 5.22, [29], for detailed descriptions, including, e.g., 
explicit computations for t^± , t^± . 

It will be convenient in constructing the Green kernel to define the adjoint normal 
modes. Thus, let us denote 

(§- $-) := $-)"^e-^ (2.19) 

We then obtain the following estimates. 

Lemma 2.5. For \p\ sufficiently small and \xi\ sufficiently large. 



i>fix,, A,0 = 72i,vi±(A,0e-'^^ (^'«)-y/(A,0(l + O(e-^l-^l)), 



(2.20) 



where fj,^ are defined as in Lemma \2.^ and Vp are dual eigenvectors of A± . Here, 
as in Lemma \2.3[. the factors 

721,^±.721,0± ~ 1 

for fast and intermediate- slow modes and for hyperbolic super-slow modes, and 

72i,^± ~ 1 + r'l^SmA - r/f (01 + p]^^ 

l2wf ~ 1 + [P~'\'^rnX - r/f (01 + p\^f 
for glancing super-slow modes, for some t^± , t^± < 1 depending on Sj . 
Proof. The proof is straightforward from the estimates of -ipf , (pf in f l2.18p . □ 
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2.3. Solutions near xi ~ 

Our goal now is to analyze system (12 ■4p close to the singularity xi = 0. To fix 
ideas, let us again stick to the case xi > 0, the case xi < being equivalent. We 
introduce a "stretched" variable as follows: 

dz 



6 



ap{z) 

so that = 0, and — )■ +oo as Xi 0+ (note that thanks to (ll.lOp . ap{z) ~ z 
when z is small). Under this change of variables we get 

, du 1 du 1 . 

dxi ap{xi) dC,i ap(xi) 

and denoting ' = d/d^i. In the stretched variables, making some further changes 
of variables if necessary, the system (12. 6 p becomes a block-diagonalized system at 
leading order of the form 



-a 




Z + ap(6)B(6)^, (2.21) 



where B(^i) is some bounded matrix and a is the (p, p) entry of the matrix \ + iA^-\- 
LB + {T~^)x^AT + Ax^, noting that due to the positive diffusion assumption (S2) on 
LB and definitions of L = T~^L and B = BT, we have 

3fJe a(6) >So>0, 

for some Sq and any sufficiently large or Xi sufficiently near zero. 

The blocks —al and are clearly spectrally separated and the error is of order 
0{\ap{^i)\) — )■ as ^1 — > +00. By the standard pointwise reduction lemma (see, for 
example. Proposition B.l, [21j), we can separate the fiow into slow and fast coor- 
dinates. Indeed, after proper transformations we separate the fiows on the reduced 
manifolds of form 

Zi = -aZi + 0{ap)Zi, (2.22) 
Z2 = 0{ap)Z2. (2.23) 

Since —^ea < —6q < for p ~ and > 1/e, with e > sufficiently small, and 
since cip(^i) — )■ as ^1 — )■ +00, the Zi mode decay to zero as ^1 — )■ +00, in view of 

These fast decaying modes correspond to fast decaying to zero solutions when 
Xi — )■ O"*" in the original w-variable. The Z2 modes comprise slow dynamics of the flow 
as Xi — )■ O"*". We summarize these into the following proposition. 
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Proposition 2.6. ^2], Proposition 2.4] There exists < eo ^ 1 sufficiently small, 
such that, in the small frequency regime A ~ 0, the solutions to the spectral system 
(I2.4p in (—60, 0) U (0, eo) are spanned by fast modes 



A) = ( gi I ± 60 ^ a:i ^ 0, (2.24) 



decaying to zero as xi — )■ 0^, and slowly varying modes 

'if 



4i:^u A) = I 5* I , ±e„ § 1-1 % 0, (2.25) 



with bounded limits as Xi — )■ 0^. 

Moreover, the fast modes fl2.24p decay as 



<.~kir ^0 (2.26) 



and ^ 

q£ j ~ Oi\x^rap{x^)) ^ 0, P, (2-27) 

as Xi — )■ 0^; here, is some positive constant, kp = n — p + 2, and 

2.4. Two Evans functions 

Having constructed bases of the solutions in regions x < —Xq, x > Xq, and |x| < 
|xo|, we can extend the modes (p^ in $^ to regions of negative /positive values of y 
by expressing them as linear combinations of solution bases constructed in Lemmas 
12.31 and 12.61 in these respective regions. Thus, we are able to define the following two 
variable-dependent Evans functions 

D+(yi, A,0 := det($+ WjT^ ^~){y,,\,0, for y, > 0, (2.28) 

and 

D_(yi, A,|) := det($+ iy+ $-)(yi, A,|), for y^ < 0, (2.29) 

where are defined as in fl2Tfll . and W^^ = {u^^, qu^.p^Y as in ^TIM . and 

kp = n — p + 2. 

We observe the following simple properties of D±. 
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Lemma 2.7. For X sufficiently small, we have 

D±iy,,X,0 = l±iyi)idetA,)-'A{X,0 + Oip'), (2.30) 
where A(A,^) is the Lopatinski determinant, defined as 

A(A,0:=det(r2+ ■■■ r+_, r^^^, ■■■ r-^, X[u] + z[f^{u)] 
7±(yi) := det 



qt ql\ (2-31) 



with [u] = u+ — U- and constant eigenvectors of {A]_) ^(LB)±, spanning the 
stable/unstable subspaces at ±oo, respectively. 

Proof. The computation follows straightforwardly from lines of the computations in 



29l . pp. 59-61], and those in Lemma 2.5]. □ 



Lemma 2.8. Defining the Evans functions 

^±(A,0:=^±(±1,A,0, (2.32) 

we then have 

D+{X,0 = mD4X,0 + 0(p') (2.33) 
for some nonzero factor m. 

Proof. Proposition 12.61 gives 



= 1 ^t, I = Oi\xn, (2.34) 



as Xi — )■ 0, where ao is defined as in Proposition l2.6[ Thus, there are positive constants 
ei, €2 near zero such that 



Thus, this together with the fact that Wj^^ are solutions of the ODE (12.71) yields 

for some nonzero constant m^^. Putting these estimates into fl2.30p and using conti- 
nuity of D± in (A,^) near zero, we easily obtain the conclusion. □ 
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3. Resolvent kernel bounds in low— frequency regions 

In this section, we shall derive pointwise bounds on the resolvent kernel Q-,^ ^{xi, yi) 
in low-frequency regimes, that is, |(A, f)| — )■ 0, following closely the analysis developed 
for the one-dimensional stability in |2l|. We recall that the Green kernel Q^^{xi,yi) 
of the ODE system f l2.7p solves 



dxr{<S>{xi)g^^^{xi,yi)) - A{xi,X,OGx,i{xi,yi) = Sy,{xi), (3.1) 

in the distributional sense, where 5yi(-) denotes the standard Dirac function with 
mass at Xi = yi. The kernel Qxii-^V^) then satisfies the jump conditions at Xi = yi. 

/A^{y,)-' 0\ 

[M-'^i)]= (3-2) 
V ij 

For definiteness, throughout this section, we consider only the case yi < 0. The case 
2/1 > is completely analogous by symmetry. 

By solving the jump conditions at Xi = yi, one observes that the Green kernel 
Q^^^{xi,yi) can be expressed in terms of decaying solutions at ±oo as follows 

^ ( h+ix,,x,Oc^iyiA,0 + w^^ix^,x,Ocliyi,x,i), ^i>yi, .33. 

''^^""''^'^ [-<^-{x,,X,OC-{yi,X,0, x,<y,. 

where are row vectors. We compute the coefficients by means of the trans- 
mission conditions (13. 2p at yi. Therefore, solving by Cramer's rule the system 

%(l/i)-i 0^ 

($+ ^ )\Cl\ =1 1 I , (3.4) 

we readily obtain. 








Cljiyi,KO = D-{yuX,0-'{^+ W^^ <^-y'' \^ 1 0| (3.5) 

where M"''^^ denotes the adjugate matrix of a matrix M. Note that 
Cf^iy,A,i)=aMr'D4y„X,i)-'{^^ (i/i, A, 0, (3.6) 

Cf,iyuX,i) = J2D-iyi,>^,ir'{^^ K ^"f {yiA,OiA,{y,r')M, l^P, 

k 

(3.7) 
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where O'-' is the determinant of the (i, j) minor, and {Ai{yi)~^)ki, I 7^ p, are bounded 
in yi. Thus, Cj^ are only coefficients that are possibly singular as yi near zero because 
of singularity in the p^^ column of the jump- condition matrix (13. 2p . 

Lemma 3.1. Define p := |(A,^)|. For p sufficiently small and for yi near zero, we 
have 

( Oip-'), j = l,n + 2, 

Ct{yi,\0 = \ 0(ap(i/i)-i|yi|--«), j = kj„ (3.8) 
y 0{1), otherwise, 

where kp = n — p + 2, ao is defined as in Proposition \2.6\ and 0{1) is a uniformly 
bounded function, probably depending on yi,X,^. 

Proof. We shall first estimate C~^2,piyi^ ■^^0- Observe that 

($+ $-r^'(l/i,A,0 = ($+ $-)""^^i/,0) + O(p) 

where by the same way as done in Lemma 12.71 we obtain an estimate 

($+ $-)^'""''(yi,0) = a,(detAi)-S-(yi)A^'"+^(A,O, 

where 7-(?/i) and A(A,^) are defined as in (I2.3ip . and AP'"+^(A,^) denotes the minor 
determinant. Thus, recalling (I2.30p and (13. 6p . we can estimate Cn+2,piyi^ ^^0 

C-^2,MA,i) = aMr'D4y,,X,0-'{'^^ W^^ <f-)''"^' (l/i, A, 
= ^AP'"+2 + (9(1) 

where 0(1) is uniformly bounded since normal modes are all bounded as yi is 
near zero. Similar computations can be done for C~_^_2 i- Thus, we obtain the bound 
for C~^2 claimed, by our strong Evans function assumption (D). The bound for 
Ci follows similarly, noting that (f)n+2 = (pt sd. p = 0. 

For the estimate on C^, we observe that by view of the definition (I2.30p of 

D_{yi, A, ^) and the estimate M/^(?/i) ~ {0, ■ ■ ■ ,\yi\'^ , ■ ■ ■ , 0)* in Proposition 12. 6[ 

\D4y,,X,i)\^\y^r\D4\,i)\. (3.9) 

This together with the fact that (f)n+2 = at p = yields the estimate for as 
claimed. 

Finally, we derive estimates for (resp. C^^) for 1 < j < hp (resp. kp < 
j < n + 2). By again applying the estimate (12.240 on Wkp and using the fact that 
0~^2 = at p = 0, we obtain 

($+ iy+ <I>~Y' = Oip\y^rapiy^)) 
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and for k ^ p, 

These estimates together with (13.91) and fl3.7l) . fl3.6p yield estimates for as claimed. 

□ 

Proposition 3.2 (Resolvent kernel bounds as — )■ 0). Let W = {U, Q) he the shock 
profile. Then, for p sufficiently small and for yi near zero, there hold 

g^;^{x,,y{) = 0{p-^){W'{xr) + pO(e-^l^-^l)) + 0{e-p'\^-\) (3.10) 

for ?/i < < xi, and 

g,,^{x^,y^) = 0(p-i)(iy'(xi) +p0(e-''l-l)) + 0(l)(l + ^^^J^) (3.11) 

for yi < xi < 0, and 

g^^^{x,,y,) = 0{p-'){W'{x,)+pO{e-''\^^\)) 
for Xi < yi < 0, for some t] > 0. Similar bounds can be obtained for the case yi > 0. 

Proof. For the case < < xi, by using the facts that M/^(x) = and that 

(l)f{xi, A, ^) is the fast-decaying mode at x = +oo which we can choose (pf{xi, 0, 0) = 
W'{xi), together with the estimate (13.81) . (13.31) then becomes 

kp-l 

g,^^{x^, yi) = $+(xi, A, OC^iyi, A, 1) = ^t(^^^ ^' OC^yi, A, 

kp — 1 

= 0{p-') (W, + O(p)e-^l^l) + 0{1) J2 

i=2 

yielding (13.101) . noticing that ^ep^ < —9p^ for j = 2, kp — 1. Similarly, the last 

case xi < 1/1 < is obtained by the estimate (13.81) and the fact that W^{xi, A,^) is 
the fast-decaying mode at x = — oo, and iy3"(xi,0,0) = W'{xi). In the second case 
yi < Xi < 0, the formula (13. 3p reads 

g,^^{x^,y,) = ^+ix^A,i)C^{ylA,i) + W+ix^,\,i)Cliy^,\,i) 

where the first term contributes the terms as in the first case, and the second term is 
estimated by and i^M)- □ 
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Next, we estimate the kernel Q,^ ^{xi,yi) for yi away from zero. We then obtain 
the following representation for ^(xi, yi), for yi large. 

Proposition 3.3. Under the assumptions of Theorem for \p\ sufficiently small 
and 1 2/1 1 sufficiently large, we have 

M^i'^i) = (3.12) 

for yi < < xi, and 
Gx,i{xuyi) = ^c/J(A,O07(a;i, A,f)V^fc (yi, A,0* - ^ ^/'^(xi, A, OV^fc (l/i, A, 1)*, 

j,k k 

(3.13) 

for yi < xi < 0, and 

Gx,i{xuyi) = ^dji^i^, {xi, X, {yi, X, ^)* + ^^^(xi, A,|)0^(?/i, A,|)*, 

j,k k 

(3-14) 

for Xi < yi <0, where c^^(A, ^), (i^^(A, ,^) are scalar meromorphic functions satisfying 

c^ = (-4, 0)(<f+ 

and 

Proof. By using the representation fl3.3p and expressing the normal modes in terms 
of the solutions in basis in each region yi > or yi < 0, the proof follows easily by 
direct computations. □ 

We define 

r«:={A : ^eX = -9i{\i\^ +\QmX\^)}, (3.15) 

for > and |(^, A)| sufficiently small. Applying Proposition 13.31 and Lemmas 12.31 
and 12.51 we obtain the following proposition. 

Proposition 3.4 (Resolvent kernel bounds as \yi\ — > oo). For A G and p : = 

|(^, A)|, 9i sufficiently small, for \yi\ large enough, there holds 

KGxM^yi)\ <C72P^(p-'e-^l^^le-^'''l^^l +e-^^'l^-^^l), (3.16) 
for /3 = 0,1, and 72 defined as 

72(A,0 := 1 + $^[p-^|53mA - r/f (01 + p]'/'^-', (3.17) 

r]f,Sj defined as in (H2). 
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Proof. The estimate fl3.16p is a direct consequence of the representation of Q;^ ^{xi, yi) 
recalled in Proposition 13.31 and the estimates on the normal modes obtained in Lem- 
mas 12.31 and 12.51 recalling the uniform Evans function condition gives |-D±|~^ = 
O(p-i). □ 



Corollary 3.5. For A G and p := |(^, A)|, 6i sufficiently small, there holds 





Xi 


V 




i)\yi 





(3.18) 



for (3 = 0,1, where x = 1 for —l<yi<xi<OorO<xi<yi<l and x = 
otherwise, and 72 is defined as in fl3.17p . 

Remark 3.6. The last term in fl3.18p accounts for the singularity of the Green kernel 
when yi is near the singular point yi = 0. 

4. Solution operator estimates 

The solution operator S{t) := e^* of the linearized equations may be decomposed 



into low frequency and high frequency parts as S{t) = Si{t) + S2{t) as in [28[, where 



|5|<r Jr« 

and 

9i+ioo 



where we recall that 

:= {A : 3fJeA = -^idlP + I^SmAH}, (4.3) 

for ^1 > sufficiently small. 

Then, we obtain the following proposition. 

Proposition 4.1. The solution operator S{t) = e"^* of the linearized equations may 
be decomposed into low frequency and high frequency parts as S{t) = Si{t) + S2{t) 
satisfying 



\SimidlfU. <C(1 +t)-^^^"VP)-T|;|^, 

-^(i-i/p) 



. w'l (4-4) 

+ C{1 + t)-^^^-^/^^-^-^|/|il(,;^,l + ..(.,)) 
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for all 2 < p < CO, d > 2, and (3 = with f3i = 0,1, where \f\Li(^x;H^+i^i{xi)) 

denotes the standard space in x and the H^^^^ Sobolev space in xi, and 



\&Ild2s,{t)f\,.<ce-'^y\^,,,,,,,, 

for 7 = (71, 7) with 71 = 0, 1. 

The following subsections are devoted to the proof of this proposition. 

4.1. Low-frequency estimates 

Bounds on Si are based on the following resolvent estimates. 



(4.5) 



Proposition 4.2 (Low-frequency bounds). Under the hypotheses of Theorem M . 1\ for 
A G (defined as in f l3.15p ) and p := |(^, A)|, 9i sufficiently small, there holds the 
resolvent bound 



|(£|-A)-^9j/U.(.,)<Cp-^+^2|/|LV.,) + CK/U-(^i), 
for all 2 < p < 00, /3 = 0, 1, and 72 defined as in f l3.17p . 
Proof. From the resolvent bound fl3.18p . we obtain 



(4.6) 



|(L^--A)-^9fjU.(.,)= / g^^,ixuyi)d;jiy,)dyi 



LP{xi) 

+ e-''"^^^-y^\)\f{yi)\dy 



LP{xi) 





Xi 


V 






V 



The first term in the first integral is estimated as 



Kf(y^)\dyi 



LP{0,1) 



72 



^'-1 p-'^'e-'p'\y^\\f{yi)\dyi 



< C72p-'+^|/|LH.i 



LP(xi) 



(^1) 



and, by using the convolution inequality \g * h\Lv < Is'IlpI^IlI; the second term is 
bounded by 

C^2P^\e-"'"\\pix,)\fWix,) < C72p-'/^+^|/|Li(.,)- 
Finally, for the last term, we use the fact that a{yi) ~ ?/i as ?/i — and j^^ {l + 
j dyi < +00, for xi G (0, 1). The estimate (14. 6 p is thus obtained as claimed. □ 
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4 ■ 2. Proof of bounds for Si (t) . 

The proof will follow in a same way as done in 28|]. We shall give a sketch here. 
Let u{xi,^,X) denote the solution of (L^ — A)^ = /, where /(xi,^) denotes Fourier 
transform of /, and 



uix,t) ■.= Siit)f 



(27ri) 



|g|<r ^r«n{|A|<r} 



Recalling the resolvent estimates in Proposition \A.2\ we have 

\u{Xi,^, X)\lp{x,) < Cp''^-f2\f\LHxi)+C\f\Hi{xi) 

< C'p~S2|/|Li(aO + C\f\L^ii.m{xi))- 

Therefore, using Parseval's identity, Fubini's theorem, and the triangle inequality, we 
may estimate 



I 1 2 / \ 

mL^xux)W = (^2Tryd 
1 

1 

< 



Xl 



e^^u{xi,^,X)dX 



r«n{|A|<r} 



2 _ 



e^^u{xi,^,X)dX 



(27r)2'i 
'li(x) 



r«n{lAl<r} 



r«n{|A|<r} 



L2(xi) 



d^dxi 
di 



JkeXt 



u{^i,iA)\L^{x^)dX 



2 _ 
d^ 



<C\f\' 



rfn{|A|<r} 



2 _ 

d^ 



+ C\f\ 



r«n{|A|<r} 



2 _ 



Specifically, parametrizing by 

xiik) = tk-eiie + \^\^), ke 

and observing that by (13.171) . 



12P 



-<(w + ifi)-[i + E(^^)""' 

j 

<(w + ki)-'[i+E(^^)""' 



P 

where e := — - — with recalling that are defined in (H2), we estimate 
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(4.7) 



'r«n{lA|<r} 



2 _ 



< le 



2 _ 



-2ei|C|2t|^|_2e 



2 _ 



,-2ei|5|2t|||-2e 



2 _ 



e-'''''\k\'-^dk 



2 _ 
di 



and 



r«n{|A|<r} 



2 _ 



-ei(fc2+|5-|2)t^^ 



2 _ 



Similar estimates can be obtained for the L°° bounds and thus the bounds 
by the standard interpolation between LP' and Also, the Xi-derivative bounds 
follow similarly by using the resolvent bounds in Proposition 14.21 with /3i = 1. The 

x-derivative bounds are straightforward by the fact that d^f = (if)^/. 

4 ■ 3. Proof of hounds for S2 (t) . 

The bounds for S2{t) are direct consequences of the following resolvent bounds. 

Proposition 4.3 (High-frequency bounds). For some R, C sufficiently large and 9 > 
sufficiently small, 

\i^^->^r'f\HH^,)<C\f\HH^,), (4.8) 

and 



C 



(4.9) 



for all |(^, A)| > R and MA > —6, where f is the Fourier transform of f in variable 
X and \f\Hi(^,) := 1(1 + \^x^\ + l^Dfli^x^)- 

Proof. The proof is straightforward by deriving an energy estimate as a Laplace- 
Fourier transformed version with respect to variables (A, x) of the nonlinear damping 
energy estimate, presented in the next section (see, for example, an analog proof 
carried out in [ij]. Section 6, to treat the one- dimensional problem). □ 



We also have the following: 
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Proposition 4.4 (Mid- frequency bounds). Strong spectral stability (D) yields 

\{C-^-\Y^\hi^^^)<C, for R-^ <\{i,X)\ < R andmX> -e, (4.10) 

for any R and C = C{R) sufficiently large and 9 = 9{R) > sufficiently small, where 
Ifln'^ixi) '^^ defined as in Proposition \4.3\ 

Proof. This is due to compactness of tlie set of frequencies under consideration to- 
gether with the fact that the resolvent (A — C^)~^ is analytic with respect to in 

(I, A). □ 

Proof of bounds for S2 (t) . The proof starts with the following resolvent identity, using 
analyticity on the resolvent set p{C^) of the resolvent (A — for all / G V{C^), 



(A - c^r'f = X'\x - c^r'qf + x-'f. 

Using this identity and fl4.2p . we estimate 

1 



S2{t)f 



{2ni) 



9i+ioo 



+ 



1 



^|CP + |9'mA|2>6»i 

X e*«-*+^*A-i(A - C^)-^C^fXxi,Od^dX 

^|^|2 + |3'mA|2>ei 



X e'^-^+^'X-^f{xi,Od^dX 



where, by Plancherel's identity and Propositions 14.31 and 14.4^ we have 

-8l+ioD 



\Si\ 



< c 



< Ce 



— 61— ioo 
-9it 



|Ar^|e^*||(A-£^-)-i£^-/\.(^;,^)MA| 



IAI-3/2 



il —lOO 



(l + |f|)lVl//i(.0 



\dX\ 



and 



'S'2|l2 < 



+ 



PV 
1 



—t/i —too 



A-^e^*rfA / e'^-^f{xi,Od^ 



Ll 



{2'kY 
<Ce-^^*l/|L?, 



PV. 



il — ir 



X-^e^*dX I e'^<f{xi,Od^ 



Ll 



(4.11) 



(4.12) 



(4.13) 
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by direct computations, noting that the integral in A in the first term is identi- 
cally zero. This completes the proof of the first inequality stated in the proposition. 
Derivative bounds follow similarly. □ 



5. Nonlinear damping estimate 

In this section, we establish an auxiliary damping energy estimate. We consider 
the nonlinear perturbation equations for variables {u, q) 



Ut 



+ ^ Aj{x)u^^ + Ldivg = - ^ Mj{x)U^^, 



(5.1) 



-Vdivg + q + V{B{x)u) = 0, 
where we have denoted 

Aj{x, t) := dfj{U + u), Mj{x, t) = dfj{U + u) - dfj{U), 

and ^ 

B{x,t) := / dg{U{xi) + su{x,t)) ds. 
Jo 

Here, the functions Aj{x, t) and B{x, t) should not be confused with Aj{xi) and B{xi) 
that used in the previous sections. The former notation is only used in this Section. 

Proposition 5.1. Under the assumptions of Theorem so long as ||m||vi^2,oo re- 
mains smaller than a small constant ( and the amplitude It/xJ is sufficiently small, 
there holds 

\u\l,{t) < e-^^'lull^iO) + C [ e-^^'-^Vli^ls) ds, 7] > 0, (5.2) 

Jo 

for k = 1, ...,s, with s large as in Theorem \l.l\ 

Proof. We symmetrize the hyperbolic system in (15.11) as 

AoUt + ^ Aj{x)u^^ + AoLdivg = - ^ Mj{x)U^^ 
i j 

where Aq is the symmetrizer matrix and Aj = A^Aj, Mj = AqMj. We then observe 
that 

\Ao,\,\Aot\A^jx\A^jt\AMj,\,\Mjt\ABx\ABt\ = 0{\U,A+C)- (5.4) 

Taking the inner product of q against the second equation in (15. ip and applying the 
integration by parts, we easily obtain 

iVgli^ + \q\h = {Bu, Vg) < i|Vg|i, + C\u\l,. 
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Likewise, we can also get for k > 1 

IqIh'' ^ C\u\ffk-i, (5.5) 

for some universal constant C. 

Taking the inner product of u against the system (15.31) and integrating by parts, 
we get 

2^(^oM,m) = -l{^i^u,u) - {U^,Mj,u) - {AoLdiYq,u) 
which together with (15. 5 p and the Holder inequality gives 

j^\u\h < C\u\h. (5.6) 
Now, to obtain the estimates (15. 2p in the case of A; = 1, we compute 

= -{{AoAjU^^ + AoLdivq)^^,Uxi^) + (0(|f4il +Oux^,Ux^) 
where, noting that AqAj is symmetric, we have 

-{AoAjU^.^^,u^^) = ^{{AoA).^M^^,u^^) = (0(|f/xj +Ou^„u^,), 

and 

-{{AoLdivq)^^,u^^) = -(AoL(divg)^^, - ((AoL)^^divg, m^J 

= -{AoLBu^^,u^^) + (0(|f/^J +C)Mxfe,Mx,) + llgll^i 

= -((AoL5)±M.„n,.,) + {Oi\U,,\+C)u,,,u,,) + Oil)\\u\\l,. 

Thus, we obtain the following first-order "Friedrichs-type" estimate 

= -{iAoLB)±u^^,u^J + {0{\U^,\ + C)u^^,u^^) + 0{l)\\u\\l2. 

(5.8) 

We quickly observe that since LB is not (strongly) positive definite, the first term 
on the right-hand side of (15.80 does not provide a full control on the norm of u. 
We shall then need to apply a so-called Kawashima-type estimate. Let us first recall 
the following well-known result of Shizuta and Kawashima, asserting that hyperbolic 
effects can compensate for degenerate diffusion LB, as revealed by the existence of a 
compensating matrix K. 



29 



Lemma 5.2 (Shizuta-Kawashima; Assuming (Al), condition (A2) is equiva- 

lent to the following: 

(Kl) There exist smooth skew- symmetric '^compensating matrices" K{^), homo- 
geneous degree one in ^, such that 

^e{A,LB\i\' - J^eA-)^ >0\e>Q (5.9) 

3 

for all^e R'' \{0}. 

We now use this lemma to give sufficient (or rather, iJ^) bounds. Let K{^) 
be the skew-symmetry from the Lemma 15.21 We then compute 

X 11 

2^(^(^^'J"'") = {Kut,u) + -{Ktu,u) - -{K^^u,ut) 

= -{KA,u,^+KLdiyq,u) + {0{\U^,\ + Qu.^u,,) + 0{l)\\u\\l2 
= -{KA,u^^,u) + (0(|f/,J +C)n,„M,,) + Oil)\\u\\l,. 
= -((irA,)±«.,,«) + (0(|t/.,| +CK„w,.,) + o(i)||«||i.. 

Using Plancherel's identity, we then obtain 

^ '^■{K{d,)u,u) = {(J2k{0^,Aj)±u,u) + {O{\U,,\+0u,,u,) + O{l) 



2dt 



(5.10) 

where u is the Fourier transform of u in x; here, stands for d^f, for some Xk- 

Let us now combine the above estimate with the Friedrichs-type estimate. By 
adding up (15. 8p and fl5.10p together, we obtain 

j 

which, together with (15. 9p and the fact that C(|f/a;| + () is sufficiently small, yields 

--((ir(9,)M,M) + {Aou,,u,)^ < --9{u,,u,) + 0{l)\\u\\l,. (5.11) 

Very similarly, we also obtain the following estimate for higher derivatives 9", |a| = 
A; > 1, 

^±(^{Kid^)d:-\ dr'u) + {Aody, 9») < -^eidy, dy) + o(i) 

(5.12) 
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To conclude the desired estimates from the above Kawashima and Friedrichs- 
type estimates, we define 

s 

fc=0 \a\=k 

for (5 > 0. By applying the standard Cauchy's inequality on {K{dx)d^~^u, d^~^u) and 
using the positive definiteness of Aq, we observe that S{t) ~ We then use the 

above estimates (15. lip and f l5.12p . and take 6 sufficiently small to derive 

±S{t)<-e,S{t) + C\\u\\Ut) (5.13) 

for some ^3 > 0, from which (15. 2p follows by the standard Gronwall's inequality. The 
proof of Proposition 15.11 is then complete. □ 

6. Nonlinear analysis 

Defining the perturbation variable u := u — U , we obtain the nonlinear perturba- 
tion equations 

ut- Cu = ^N^{u,Ucc)xj, (6.1) 
j 

where N^{u,Ux) = 0{\u\\ux\ + \u\'^) so long as |m| remains bounded. We then apply 
the Duhamel formula (I1.16P to (16. ip . yielding 



u[x. 



t)=S{t)uo+ [ S{t-s)y^dx^N^{u,Ux)ds (6.2) 

"'0 



where u{x,0) = Uo{x), recalling that S{t) = e^^ denotes the linearized solution oper- 
ator. 

Proof of Theorem Define 

C(t) ■= sup (\u{s)\l2{1 + s)^ + \u{s)\l^{1 + s)"^-') (6.3) 

where e > is arbitrary small in case of c? = 2 and e = in case of (i > 3. 

We first show that C,{t) is well-defined at least locally in time. Indeed, the sym- 
metrizability assumption (Al) easily yields the following a priori H'^ "Friedrichs-type" 
estimate (see also (15. 8 p for an version): 



d 



umis < c\\u{t)\\is[i + \Ht)\u^), 



dt 
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for some positive constant C and s > 1 + d/2. It is then easy to see that the 
standard short-time theory and local well-posedness in if* can be applied for the 
perturbation equations f l6.ip . from a standard nonlinear iteration scheme and the 



above a priori estimate. See, for example, [30|, Proposition 1.6, for a detailed proof 
of the local well-posedness for symmetrizable hyperbolic and hyperbolic-parabolic 
systems. Furthermore, the local-wellposedness argument also shows that the solution 
u G H'^ indeed exists on the open time-interval for which \u\h'> remains bounded, and 
thus on this interval C,{t) is well-defined and continuous. 

We shall prove next that, for all t > for which the solution exists with (^{t) 
uniformly bounded by some fixed and sufficiently small constant, there holds 

m<C{\u,WnHs+atf). (6.4) 

This bound together with continuity of C(t) implies that 

at) < 2C|Mo|LinH= (6.5) 

for t > 0, provided that iMolLin/f^ < 1/4C^, by the standard continuous induction 
argument. Indeed, assume that (16. 5p fails. By continuity, we can take the first T > 
such that C(^) = 2C|Mo|LinH=- The estimate (16. 4p then yields 

= C\uo\LinH- (l + 4C^|Mo|LinH=) • 

A contradiction then occurs if the initial perturbation is small, namely IuoIlidh" < 
1/4C2. 

In addition, we observe that the claim also provides sufficient bounds on if* norm 
of the solution. To see this, we apply the Proposition 15.11 and the Sobolev embeding 
inequality |m|vi^2,cx) < C\u\h'>- We then have 

\u{t)\ls < Ce~''\uo\jjs + C f e-'^'-^^uiT)\l2dT 

Jo 

< c{\uo\is+aty){i+tp'-'y\ 

With such a uniform bound on if* norm, the solution can then be extended to a 
larger time interval. Repetition of these arguments yields the global existence of the 
solution, provided that the claim (16. 4 p is proved uniformly in time. This and the 
estimate (16. 5p would then complete the proof of the main theorem. 
Thus, it remains to prove the claim (16. 4p . First by ()6.2p . we obtain 

\u{t)\L2 < \S{t)uo\L2+ [ \Si{t - s)d,^Nas)\L2ds + [ \S2{t-s)d,^N^{s)\L2ds 



(6.6) 







h + h + h 

(6.7) 
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where by using the estimates in Proposition 14.11 we estimate 
Ji := \S(t)uo\L2 < C(l + t)"^|uo|LinH3, 
\Si{t~s)d,^N^{s)\L2ds 

< C [ {l + t-sr"-^--^{\N^{s)\Li + \d^,N^{s)\Lri,-,HH.,)))ds 
Jo 

< C {1 + t- sy'^'^luljjsds 
Jo 

< ci\uo\is+atf) [ + 

Jo 



and 



< c{i+tr^{\uo\is+atr) 



= [ \S2{t-s)d^^N^{s)\L^ds < [ e-'^'-'^\d,^N^{s)\Hsds 
Jo Jo 

< cf e^^^'~'\\u\Lo^ + \u^\L^)\u\H^ds < c[ e~^^'-'^\u\lsds 

Jo Jo 



< C{\uo\ls+atY) / e-'^'-'\l + s)-^ds 

Jo 

< c{i+t)-'^{\uo\%+cm- 



Combining the above estimates immediately yields 

+ t)^ < CHuoIlw + CW')- (6.8) 

Similarly, we can obtain estimates for \u{t)\L^, noting that a Moser-type inequality 
(precisely, Lemma 1.5 in js^]) is used to give: \N{t)\i^ < C\u{t)\'jjs. This then 
completes the proof of the claim (16. 4p . and therefore the main theorem. □ 
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